function y = funcKxR2(nx, nz, zeta, CONSTS)

    [~ KxS2] = funcKxS(nx, nz, zeta, CONSTS);

    Kx = funcKxfull(nx, nz, zeta, CONSTS);
    
    y = Kx - KxS2;

end